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Abstract 

We develop a McKean-Vlasov interpretation of Navier-Stokes equations with exter- 
nal force field in the whole space, by associating with local mild L?"— solutions of the 
3d- vortex equation a generalized nonlinear diffusion with random space-time birth that 
probabilistically describes creation of rotation in the fluid due to non-conservativeness 
of the force. We establish local a well-posedness result for this process and a stochas- 
tic representation formula for the vorticity in terms of a vector-weighted version of its 
law after its birth instant. Then, we introduce a stochastic system of 3d vortices with 
mollified interaction and random space-time births, and prove the propagation of chaos 
property, with the nonlinear process as limit, at an explicit pathwise convergence rate. 
Convergence rates for stochastic approximation schemes of the velocity and the vorticity 
fields are also obtained. We thus extend and refine previous results on the probabilistic 
interpretation and stochastic approximation methods for the non-forced equation, gen- 
eralizing also a recently introduced random space-time-birth particle method for the 2d 
Navier-Stokes equation with force. 
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1 Introduction 

The Navier-Stokes equation for an homogeneous and incompressible fluid in the whole plane 
or space, subject to an external force field F, is given by 

1^ (u • V)u = z^Au - Vp F; 

(1) 

div \i{t,x) = 0; u(t,x) — > as |x| — > oo. 

Here, u denotes the velocity field, p is the (unknown) pressure function and > is the 
(constant) viscosity coefficient. When F = or, more generally, when F = V^' is a conser- 
vative field, a probabilistic interpretation of equation ([1]) in space dimension two was first 
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developed in 1982 by Marchioro and Pulvirenti [19]. Their approach was based on the vortex 
equation satisfied by the (scalar) field curl u, which in 2d and for the case of a conservative 
external field, was interpreted as a nonlinear Fokker-Planck (or McKean-Vlasov) equation 
with signed initial condition. This was associated with a nonlinear diffusion process in the 
sense of McKean, involving singular interactions through the kernel of Biot-Savart. (For 
general background on the McKean-Vlasov model, we refer the reader to Sznitman [25] and 
Meleard [20]). This approach led them to the definition of a stochastic system of particle 
or vortices with "mollified" mean field interaction, for which the time-marginal empirical 
measures converge to a solution of the vortex equation associated with ([T]). The conver- 
gence on the path space of that particle system (or equivalently, the propagation of chaos 
property) was proved later by Meleard in [21]. Those works provided a rigorous mathemat- 
ical meaning of Chorin's vortex algorithm, heuristically proposed in [1] as a probabilistic 
method to simulate the solution of the 2d-Navier-Stokes equation (see also [5j ) . 
In dimension 3, the vorticity field w = curlu is a solution of the vectorial nonlinear equation 

^ + (u • V)w = (w • V)u + z^Aw + g, 

(2) 

div wq = 0, 
where g = curl F and where the relation 

u{t,x)=K{y^){t,x):=-^ [ j^^l^l Aw{t,y)dy (3) 

holds, thanks to the incompressibility condition divu = and the so-called Biot-Savart law. 
Here, A stands for the vectorial product in M'^, K{x)A := —-^j^s^ is the three dimensional 
Biot-Savart kernel and K is the Biot-Savart operator in 3d. (We refer to Bertozzi and 
Majda [I] for this and for background on vorticity.) 

In absence of external forces, the problem of proving the approximation of solutions of the 
3d Navier-Stokes equations by a stochastic system of mean field interacting particles was 
first addressed by Esposito and Pulvirenti [H|. In that work, an approximation result of 
local solutions by a stochastic system of three dimensional vortices with cutoff and mollified 
interactions was obtained for each time instant, for initial vorticities that belonged to 
together with their Fourier transform. The convergence held for mollifying parameters that 
depended on the realizations of the empirical measures of the paths of the driving Brownian 
motions. 

Recently, we considered in [10] the mild version of the 3d-vortex equation in the "super- 
critical" spaces p > |, in the case g = 0. We proved local (in time) well-posedness 
and regularity results for that equation and, under an additional assumption on wq, we 
showed the equivalence between such solutions and a generalized nonlinear McKean-Vlasov 
process with values in X M38)3 

and singular drift term at i = 0. We then introduced a sys- 
tem of stochastic 3d vortices with cutoff and mollified interaction and proved the pathwise 
propagation of chaos property, with as limit the nonlinear process, and deduced stochastic 
particle approximation results for the velocity and vorticity fields. (We refer to [TT] for a 
rectification of the discussion in [10] about the work [8].) 

During the preparation of this work, we have become aware of the more recent work of 
Philipowski [22], who obtains (also in the case g = 0) a convergence rate for a mean field 
particle approximation of the vorticity field, for a simpler variation of the system introduced 
in [To]. (The propagation of chaos property was not addressed.) 
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In presence of an external force field, the additional additive term g = curlF in the (2d 
or 3d) vortex equation is physically interpreted as creation of rotation in the fluid. In 
order to describe this phenomenon probabilistically, a nonlinear McKean-Vlasov diffusion 
process with random space-time birth was recently associated with the 2d- vortex equation 
in Fontbona and Meleard [12j. More precisely, the law Po{dx, dt) of the instant and position 
of birth was suitable defined in terms of the initial vorticity and of the external field curlF, 
and it was shown that a scalar-weighted version of the time marginal law of this process 
after its birth time was equal to the solution to the 2d vortex equation (with data) in 
a given interval. The propagation of chaos property was established for an approximating 
system of interacting vortices, which were given birth independently at random positions 
and times following the law Pq, and a pathwise convergence rate was obtained under slight 
additional integrability assumptions on the data. 

The first purpose of the present paper is to extend the results of [lOj and [12j to the 3d 
Navier-Stokes equation with non conservative external force field. More precisely, let us 
assume that tt;o : K'^ ^ and g : x [0, T] are L^-fields. Denote by I3 the identity 

matrix in M.^ and let (Bt) be a standard 3d Brownian motion. Our main goal will be to 
study the well posedness of the following nonlinear process with singular interaction kernel 
and values in x M^®^, 

Xt = Xo + V2^^ /o l{s>r}dBs + J^K{^^^){s,Xs)l{s>r}ds, 

t£[0,T] (4) 

$j = /3 + /o VK(w)(s,X,)$,l{.>r}t^S- 

Here, {t,Xq) is a random variable in [0,T] x (independent of B) with law 

Po{dx,dt) oc 5o{dt)\'Wo{x)\dx + \g{t,x)\dx dt, 
w is the vectorial measure defined from the law of (r, X,^) by the relation 

E{i{Xt)<^th{T,Xo)l{t>,y)= [ ¥{y)y,{t,y)dy, f-.R^^R^ (5) 

and /i in ([5]) is the density with respect to Pq of the vectorial measure 6o{dt)wo{x)dx + 
g{t, x)dx dt. 

As we shall see, the relation ([5]) defines an equivalence between mild solutions w of (l2|) 
in the space [LP(M^)]^ n [L^(R^)]^ with p > |, and suitable strong solutions of ([i])-©. 
This correspondence extends the representation for w obtained in [10] when g = (or 
equivalently, when r = 0). More precisely, the function h : R? x [0, T] — > of ([5]) will 
describe the sense of rotation of a vortex at its random instant and position of birth (r, Xq) 
whereas, in the same way as in [10], the matrix process <I>j will account for the "vortex 
stretching" phenomenon proper to dimension 3. 

We will adapt the ideas and analytic techniques in [10] to establish local well-posedness and 
regularity results for the mild formulation of the vortex equation. Based on this, we shall 
then prove local (i.e. for small enough T > or data (ifoig)) pathwise well-posedness for 
the nonlinear stochastic differential equation (HI)-®, which will have singular drift terms 
at t = 0. 

We shall then introduce a stochastic system of n particles in R^ x R^®^ (or 3d- vortices) with 
cutoff and mollified interaction kernels, and with random space-time births. The second 
goal of this paper will be to prove the strong pathwise convergence of each of these particles 
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as n goes to oo, towards the nonlinear process, at an explicit L - rate. To that end, we 
will improve the techniques used in [lOj to study the nonlinear process (which relied on 
tightness estimates for approximating processes and martingale problem characterization). 
More precisely, by a fine use of regularity properties of the equation, inspired in ideas 
introduced in [12], we will show that the approximating "mollified processes" converge 
pathwise at the same rate at which mollified versions of the vortex equation converge to 
the original one. We will be able to exhibit that rate for a large class of mollified kernels, 
thanks to classic regularization techniques in Raviart [23] (which are also similar to those 
used in [22]). These results will imply the propagation of chaos in a strong norm and, 
classically, and an explicit rate in some pathwise Wasserstein distance W. From this we 
will also deduce convergence rates for approximation schemes of the vorticity and velocity 
fields. Unfortunately, the mollifying parameter will be required to go very slowly to as n 
goes to oo, which will yield a very slow (but not necessary optimal) rate for the particles 
convergence. 

Finally, we point out that our regularity results on the mild equation in will ensure that 
the stochastic flow 

6,t(x) = x + V2^{Bt -Bs) + u(r, 6,r(a;))dr, (6) 
is of class C^(M^), and so one can write 

{Xu ^t)l{t>r} = (er,t(Xo), V,er,t(^0))l{t>r}. (7) 

Equation ([5]) can thus be seen as a stochastic analog for the 3d Navier-Stokes equation of 
the "Lagrangian representation" of the vorticity of the 3d Euler equation v = Q (see e.g. fG] 
Ch.l), an analogy established in [8], [10] when g = 0. Lagrangian representations of the 3d- 
Navier-Stokes equations as stochastic analogues of representations for the Euler equation 
have been studied by several authors, some of which leading to (local) well posedness 
results for the equation. See e.g. Esposito et al. [7] and, for more recent developments, 
Busnello et al. [3] and Iyer [15j. The latter works follow approaches in some sense "dual" 
to ours, establishing representations of strong solutions of the vortex or Navier-Stokes 
equations in terms of expectations of the initial data, after being transported and modified 
by the stochastic flow. A related stochastic approach is adopted in Gomes [14J to establish 
a variational formulation of the Navier-Stokes equation, analog to Arnold's variational 
characterization of the Euler equation. A seemingly very different further probabilistic 
point of view, providing global well posedness for small initial data, was introduced by Le 
Jan and Sznitman in ^18j, who associated with the Fourier transform of the velocity field 
a multitype branching process or stochastic cascade. See e.g. Bhattacharya et al. [2j for 
more recent developments in that direction. 

This work is organized as follows. In Section 2 we first present a weak formulation of 
dS]) in terms of a nonlinear martingale problem, and discuss its connection with equation 
(l2|). In Section 3, we shall obtain local well-posednes (for small enough T > or data 
(zi;o,g)) and regularity results for the mild version of the vortex equation in for p £ 
(|,3) "supercritical". In Section 4 we state some results about a nonlinear Fokker-Planck 
equation with external field associated with the process with random space-time birth X 
in We use this and the previous results to show strong local- in-time well-posedness 
for the nonlinear stochastic differential equation (jl])-©. We moreover obtain the pathwise 
convergence result and estimates for approximating mollified versions of that problem. In 
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Section 5, we introduce the system of 3d stochastic vortices with random space-time birth, 
and deduce the propagation of chaos property and its rate. We also prove approximation 
results for the velocity and the vorticity of the forced 3d-Navier-Stokes equation with their 
corresponding convergence rates. In last Section 6 we shall discuss how these rates of 
convergence are slightly improved when Sobolev regularity of the initial condition and 
external field is assumed. 
Let us establish some notation: 

- By Meas^ we denote the space of measurable real valued functions on [0, T] x M'^ 

- C"*^'^ is the set of real valued functions on [0, T] x M'^ with continuous derivatives up to 
the first order in t G [0, T] and up to the second order in x G M?. C^'^ is the subspace 
of bounded functions in C^'^ with bounded derivatives. 

- P is the space of compactly supported functions on M'^ with infinitely many deriva- 
tives. 

- For all 1 < p < oo we denote by the space LP(M'^) of real valued functions on . 
By II ■ lip we denote the corresponding norm and p* stands for the Holder conjugate 
of p. We write H^^'P = pr^'P(]R^) for the Sobolev space of functions in LP with partial 
derivatives of first order in LP. 

- If is a space of real valued functions (defined on or on [0, T] x M^), then the 
notation is used for the space of M'^-valued functions with scalar components in 
E. If E has a norm, the norm in i?3 is denoted in the same way. 

- For notational simplicity, if f , g : ^ are vector fields and Z : ^ M?®^ is 
a matrix function, we will write fg := fjgi and fZ for the row-vector {i*Z)i := 
Yl^=i^j^j,i- By we denote the gradient of f, that is the matrix (Vf)ij := 

We will simply write (Vf)g for the column-vector §§'Sj)i (instead of the usual 
"(g-V)f"). 

- C and C (T) are finite positive constants that may change from line to line. 

2 The weak 3d- vortex equation and a probabilistic interpre- 
tation of the external field 

Let us recall a that vector field u) : ^ with components in V, and such that 
/r3 V f{x)w{x)dx = for all / G P, is said to have null divergence in the distribution sense. 
We write it div w = 0. 

If the following two conditions hold, we shall say that wq -.M.^ —>■ and g : M+ x — >• 

satisfy hypothesis 

(Hp): 

• There exists p G [l,oo[ such that wo G {LP{R^)f and g(t, •) G (I;'(M3))3 for all 
t G [0,r], and supjg[o,T] l|g(*, OIIp < oo- 

• div Wo = and div g(t, •) = for all t G [0, T]. 
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A necessary assumption for our probabilistic approach will be that (Hp) holds with p = 1. 
We then denote 



|g||i,T : = 



|g(s, x)\dx ds. 



In that functional setting, the following notion of solution to ([2]) will appear to be natural: 



Definition 2.1 Let wq and g satisfy (Hi). A function 

w G L°°([0, T], L^(R^))^) is a weak solution on [0,T] of the vortex equation with initial 
condition wq and external field ^ (or "weak solution") if 



i) For i,j, k = 1,2, 3, 



[0,T]xR3 



|wj(t, x)||K(w)j(t, x)|(ix dt < oo. 



[0,T]xK3 



|Wi(t,x)| 



aK(w)j 



dxi 



{t,x] 



dx dt < oo. 



(8) 



ii) For any f e {Cl^%{[0,T] 



i{t,y)w{t,y)dy = {{0,y)wo{y)dy + 



Jr^ 



f{s,y)g{s,y)dy ds 



+ 



— (s,y) + uAi{s,y) 



+Vf (s, y)K(w)(s, y) + i{s, i/)VK(w)(s, y) 



w{s,y) dy ds. 



(9) 



Remark 2.2 We observe that for function v : M"^ ^ m L^, the functions K(v) and 
VK(v) are defined a.e on W^. Indeed, the first one can be bounded by a (scalar) Riesz 
potential operator (see Stein [24^), and thus belongs to a suitable weak-L^ space. The 
second one is defined through a singular integral operator acting on v (see e.g. 11] for this 
fact), and this implies (see also [24] ) that it is an almost everywhere defined function of 
some other weak-L^ space. 

We next introduce the central probabilistic objects we shall be dealing with, which extend 
the ideas introduced in two dimensions in [12j . 

Definition 2.3 We write Ct := [0,T] x C([0,r],M3 x R^^'-^). The canonical process in Ct 
will be denoted by (r,X, $), and the .space of probability measures on Ct is written V{Ct). 
For an element F G V{Ct), we write P° = law{X) for the second marginal and P' = law{(^) 
for the third marginal. 



We shall also denote 



We then define a probability measure Po{dt, dx) on [0, T] x M'^ by 

Poidt, dx) = 5o{dt)wQ{x)dx + g(t, x)dx dt, (11) 
together with the vectorial weight function 

h{t,x) = 1{*=0}^^ (Ikolll + llglkT) + (llt^olll + ||g||l,T) l{t>0} (12) 

where 1 denotes the indicator function and the convention = 0" is made. We notice 
that |ft-(t,x)| = llwolli + lls'llijT or 0. Moreover, we have 

Remark 2.4 For measurable bounded functions f : [0, T] x — > R'^, we have 
?{t,x)h{t,x)Po{dt,dx) = / i{0,x)wo{x)dx 

[0,T]xR3 JK3 



+ / {{t,x)g{t,x)dx dt. 

J[0,T]xM3 



Consider now Q G ViCr) such that for all G [0,r], E'^{\^t\) < oo. Then, we can associate 
with Q a family of M'^-valued vector measures {Qt)t<^[o,T] defined for all bounded 

measurable function f : ^ M'^ by 

Qtif) = E'3{f{Xt)^tHT,Xo)l{r<t})- (13) 
Moreover, Qt is absolutely continuous with respect to Q^, with 

^(a;) = E'^{'^th{T,Xo)l^,^,}\Xt = x), (14) 

and its total mass is bounded by (Utt'oHi + l|g||i,T)^'^(|*^'t|)- 

Definition 2.5 We denote by VhiCr) the subset of probability measures Q S V^Ct) under 
which the process $ belongs to L°°([0,T] x Q,dt(Si P)- 

Then, we consider the following nonlinear martingale problem: 
(MP): to find P G VbiCr) such that 

• Po[t, Xq)~^ = Pq, and Pt has bi-measurable densities (t, x) p{t,x) 



fit,Xt) - f{r,Xo) - [§(s,X,) + z.A/(s,X,)ds + K(p)(s,X,)V/(s,X, 



ds. 



ds 

1 2 

< t < r, is a continuous P-martingale for all f €z C^^' w.r.t. the filtration 

• $4 = I3 + VK{p){s, Xs)'^sds, for ah < t < T, P-almost surely. 

The following statement partially explains the relation between (MP) and equation ([2]), 
and will be useful later on: 
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Lemma 2.6 Assume that the problem (MP) has a solution P S VbiCx) satisfying 



E(^j^ \^{P){t,Xt)\dt^ < oo 



and 



e(^J^ \VK{p){t,Xt)\dt^ <oo. 
Then, w := p is a weak solution of the vortex equation with external force field 



(15) 
(16) 



Proof. The assumptions on P imply that point i) in definition [2TT] is satisfied and, moreover, 
that J^'K{p){s, Xs)ds and Jl^VK{p){s,Xs)ds are both processes with integrable variation 
(and thus absolutely continuous on [0,T]). Since under P the process is almost surely 
bounded in [0,T], it follows that it has finite variation too. Thus, Ito's formula shows that 



f{t,Xt)<!>t-i{T,Xo)- f 

Jo 



ds 



is,X,) + iyAiis,Xs) + 



Vf (s, Xs)Kip){s, X,) + f (s, Xs)VKip){s, X, 



^sl{s>r} ds 



1 2 

is a local martingale for all f G C'^'a- Moreover, it follows from the assumptions on p and 
the fact that $ is bounded that it is a true martingale. Consequently, since /i(t, Xo)1|t-<(} 
is measurable with respect to .Fq, the process 



f(t,X4)$i/i(r,Xo)l|,<t} -f(r,Xo)/i(T,Xo)l|,<t| - f 

Jo 



ds 



(s,X,)+z.Af(s,X,) 



+ Vf (s, X,)Kip)is, Xs) + f (s, Xs)VK{p)is, Xs) 



$s/i(T-,Xo)l{^<4 ds 



has vanishing expectation. We take expectation and use Fubini's theorem, and we conclude 
by Remark 12.41 and the definition of p. 

□ 

The proof well-posedness of problem (MP) will be based on analytical results about the 
"mild form" of the vortex equation ([2]), which we state in next section. These will in 
particular provide a framework where the conditions required in Lemma 12.61 will hold. 

3 The mild vortex equation in with an external field 

We denote the heat kernel in by 



GUx) :=(4vr^t)-texp(-g^), 



(17) 



where v > 0. 



Lemma 3.1 For all p G [l,oo], r > p and w £ L^, there exist positive constants Co(p;r) 
and Ci{p; r) such that for all t > 0, 



i) \\Gt *u)\\r < Co{p;r)t ^'\\w\ 
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ii) WVC^ *w\\r < Ci{p;r)t 




\p- 



Proof. Use Young's inequality and the well known estimates 




t \\m 



t 2Tr, < OO. 
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□ 



Definition 3.2 Let wq and g be functions satisfying (Hp) for some p G [1, oo]. A function 
w G L°°([0, T], LP(]R^))^) is a mild solution on [0,T] of the vortex equation with initial 
condition wq and external field (or "mild solution") if 

i) The functions K.{w)i{t, x) := K-iwit, ■))i{x), i = 1,2,3 are defined a. e. on[0,T]xM.^ 
and satisfy the integrability conditions ([8]). 

ii) For dt— almost every t, the following identity holds in LP{dx) : 



We shall state in Theorems 13.51 and 13 . 71 below the analytical results we need about (jlSp . As 
we shall see, that equation will admit an abstract formulation which is the same as in the 
case g = 0, and so we will be able to adapt the techniques in [10] with no difficulties. We 
therefore provide an abbreviate account of these results. 

We simultaneously deal with a family of "mollified" versions of equation (jlSp. Consider a 
smooth function (/? : M'^ ^ M satisfying 

i) /r3 Vix)dx = 1, 

ii) /r3 \x\\ip{x)\dx < oo, 

which is called a "cutoff function of order 1". For e > 0, let : ^ M denote the regular 
approximation of the Dirac mass fe{x) = p-v3(f ). We define the convolution operators. 



where Kf, := ipe*K = K(99£). The fact that Kf, is a regular function will follow from part ii) 
in Lemma 13.31 below. To unify notation, we also write Kq = K and 'K.^{w){x) := 'K.{w){x). 
We introduce the family {B^}£>o of operators (formally) defined on functions w, v : [0, T] x 
m? ^ by 




(18) 



Wj(s,y)K(w)(s,y) dy ds 




(19) 




(20) 



[K^(w)j(s,y)v(s,y) - Vj(s, y)K^(w)(s, y)] dy ds. 
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We are interested in the following family of "abstract" equations, for e > 0: 

V = Wo + B^(v,v), 

where ^ 

wo(t,x) ■.= G'^ *wo{x)+ [ G,^_, *g(s,-)(x) ds. 

Jo 

For a given time interval [0, T] we shall work in the Banach spaces 

Fo,r,(T;p)5 Fi ,.,{T;p)) Fo,p,T and Fi^p^T, 
with norms respectively defined by 

• ll|w|||o,r,(T;p) := SUp t^^P -'\\w{t)\\r, 
0<t<T 



(21) 



\\M\l,r,{T;p) ■= sup < t^'-P r^\\^ 
0<t<T 



and 



_)_ i2^2'.p rl \ 

k=l 



dw{t) 



dxk 



• |||w|||o,p,r :- |||w|||o,p,(T;p), 

• |||w|||i,p,T := |||w|||^p^(7-.p). 

The following continuity property of the Biot-Savart kernel is crucial: 

Lemma 3.3 Let 1 < p < 3 he given and g € (|, oo) he defined \ = ^ ~ \ - 

i) For every G Lg, the integral (jl9p is absolutely convergent for almost every x and 
one has K^(ty) € L3. There exists further a positive constant Cp^q such that 



sup ||K^(i(;)||g < Cp^q\\w\\p 

e>0 



(22) 



for all w £ L-n. 



ii) If moreover w G W^'^, then we have 'K^{w) G W^''^ , with ■^'K^{w) = K y-g^ 
and 



e I dw 

k 



sup 

e>0 



dxk 



< a 



PA 



dw 



dxi 



(23) 



for all k = 1,2, 3. 



Proof. See Lemma 2.2 in for the case e = and Remark 4.3 therein for the general 
case. [] 

Lemma 3.4 i) Let p G [1, 3) and assume (Hp). Then, we have for all r G [p, ^3^) that 



Wo e -Fi {T;p) , with III Wo III 1 ^r, {T;p) < ( 1 1 "^0 I |p + ^ ||| g ||| o,p,T ) 

for some finite constant C{r,p) > 0. 
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ii) Let |<p<3, p < I < min{^^,3} and ^ <l'< Then, there exists a finite 

constant Ci{l,l';p) not depending onT > such that for all w, v G ; (^.p), 

sup|||B''(w,v)|||i^,,^(5n.p) < Ci{l,l ■,p)T 2p III will i^i^('r.p) |||v|||i^i^(r.p), 

£>0 

where 1 — ^ > 0. 

Proof. Part i) follows from Lemma |3. 11 To bound the time integral we use moreover the 
fact that for all r > p, on has 



Gr_,*g(s,-) ds 
_ 1 

On the other hand, since t^t ^ 

v( tcU^^is.-) ds 



< Ct^' 



r p 



sup US'* lip 
,te[o,r] , 



^ 2 r p> 



p ' is integrable in if and only if r < , we have 



<C't5 + |(f-F) sup ||g(t,-)llp 



from where the statement follows. Part ii) uses Lemma 13.31 and is proved in parts ii) and 
iv) of Proposition 3.1 in [lOj. See also Remarks 4.3 and 6.3 therein for the uniformity (in 
e > 0) of the bounds. [] 

Observe that the previous lemma in particular implies (taking p = r = I = I') that for 
p G (|,3), the abstract equation (pT]) makes sense in Fi^p^x for each e > 0. Now we can 
state the extension of Theorem 3.1 in [10] to the 3d vortex equation with external field: 



Theorem 3.5 Assume that (Hp) for some | < p < 3. 

a) For each T > and e > 0, equation (|2ip has at most one solution in Fo,p,r- 

b) There is a constant Tq{p) > independent of e > such that for all T > 0, wq and 
g satisfying 

T^~^ i\\wo\\p + TMo,p,e) <roip), 

each one of Equations \21\) with e > 0, has a solution G Fi^p^y. Moreover, we 
have 

sup |||w''|||i,p,T < 2|||viro|||o,p,T- 

£>0 



Proof. For later purposes, we give in details the argument of [10]. By Lemma l3.1 H i) (with 

6-p 



p in the place of r and in that of p) and Lemma 13.31 ij, we have for all v, w G Fo,p,T 



that 



It follows that if w and v are two solutions, one has 



B^(w,v)(t)||p<C / (t-s)-2F||w(s)||p||v(s)||p<is. 



|w(t) -v(t)||p < C7(|||w|||o,p,t + |||v|||o,p,t) / {t-s) ||w(s) - v(s)||p ds, 
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ds 



and iterating the latter sufficiently many times (using the identity ^{t — s) 

Ct^+^-^ for e,e>0) we get ||w(t) - v(t)||p < C /g ||w(s) - v(s)||p ds and Gronwall's lemma 
concludes the proof. 

b) We notice that for T > small enough, one has 

iC{p,p)Ciip,p;p)T^-^i\\wo\\p + TMo,p,T) < 1, 

where C{p,p) and Ci{p,p;p) are respectively the constants in parts i) and ii) of Lemma 
13.41 with all parameters equal to p. From this and Lemma 13.41 i), the same contraction 
argument used in Theorem 3.1 b) of [10] can be applied here in the space Fi^p^T- [] 

The regularity properties of the functions we shall need rely on continuity properties of 
the "derivatives" of the (mollified) Biot-Savart operator: 

Lemma 3.6 Let 1 < r < oo. 

i) For all w £ U^, we have -^'K'^{w) G Lg for k = 1,2, 3. There exists further a positive 
constant Cp depending only on r such that 



sup 

e>0 



dxi 



< Cr\\w\ 



(24) 



for all j = 1,2,3, where }i.^{w)j is the j-th component o/K( 



w 



a) If moreover w S W^'^ , then we have g^K^(t(;) E Wr 



l,r 



with 



_d_ 



K-(^^) and 



sup 

£>0 



dxidxi 



< Cr 



dw 



dxi 



(25) 



for all i,k = 1, 2, 3. 



Proof. See Lemma 3.1 and Remark 4.3 in fTOl. 



□ 



Theorem 3.7 For p G (|, 3), let S Fi^p^x, s > be the solution of (pT|) given by Theo- 
rem \3.5\ and write u^{s,x) := K^(w^)(s, x). LetC^ denote the space of Holder continuous 
functions M'^ — > of index a G (0, 1). 

i) For all r G [p, -0^), we have 

sup Ifw'^ll]^ r (y-p) < 

£>0 ' ' ' 



ii) We have 



sup sup t2 |||u^(t)||oo + ||u^(i)|| 2p-3 I 

£>ote[o,T] c p J 



2p-3 > < oo. 



(26) 



Hi) For all r G (3, ^^), i = 1, 2, 3 we have 



sup sup t2 2'^p" 
£>0 tG[0,T] 





du%t) 


+ 


du%t) 




dxi 


oo 


dxi 





< oo. 



(27) 
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In particular, the functions 



and t 



dujt) 

dxi 



i = 1,2,3 



belong to L^{[0,T], 



Proof. Observe that parts i ) and ii ) of Lemma 13.41 provide an estimate of the form 

ll|wll|i,r,(T;p) < c(/',p)(||^/;o||p + r|||g|||o,p,T) + A(r,/,r)A', 

for suitable / and and with A(T, /,/') an uniform upper bound for the norms of the 
operators B"^ : (F;^ ; c^.p))^ — > Fj^ ;/ c^.p) and Ai a given upper bound of |||w'^|||;^ ; (j^.p). Then, 
starting from the fact that the functions w"^ S Fi p (^rp;p) = Fi,p,T are uniformly bounded 
in e > 0, we can apply several times Lemma 13.41 and the previous inequality, (using also 
the fact that wq G F^^^/ (^.p) for all I' G [p, ^^)), and obtain an increasing sequence I' = In 
such that Iq = p, In (3, ^^)) and G F^ (j-.p) with liW^li ;^ (j^.p) controlled in terms 
of |||w^|||i (j'.p) and |||wo||||||i_;^^(T;p)- One can thus chose large enough such that /at > r 
and conclude with an interpolation inequality in the spaces F^ ; (2";p)- We refer to the proof 
of Theorem 3.2 ii) in [1(5] for this and for an explicit construction of the sequence In- 
Next, Lemma 13.31 and Theorem 13.51 imply that for q = -0^ > 3, 

sup lu^li^^^T < Csup Iw^lli^p^T < C"(||wo||p + r|||g|||o,p,T)- 

e>0 e>0 

Using the continuous embedding of W^'^ into n C^~m for all m > 3, we deduce part 
ii), taking m = q. To prove part Hi) we use part i), Lemma 13.61 and the same embedding 
result as before but with m = r. See Corollary 3.1 in [TO] for details. [] 



4 The nonlinear process 

Definition 4.1 VJ^s is the space of probability measures P G VbiCx) satisfying the follow- 
ing conditions: 

• For each t G [0, T], the time marginal is absolutely continuous with respect to 
Lebesgue's measure, and its family of densities p{t,x) has a version that belongs to 
Fo,p,T for some p > | • 

• div pt = for all t G [0, T] . 
In this section, we shall prove 

Theorem 4.2 Assume that (Hi) and (Hp) hold for some p G (|, 3). Then, the following 
hold: 

a) For every T > 0, the nonlinear martingale problem (MP) has at most one solution 
P in the class VTa ■ Moreover the function defined in terms of P by 

"'2 

w{t,x) := p{t,x)E^{'^th{T,Xo)l{t>T}\Xt = x) 
is then the unique solution of (jlSp in Fo,i,t H Fo,p,T- 
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b) Consider in a filtered probability space {Q,J-',J-t,P) « standard Brownian motion B 
in M^, and an measurable random variable (r, Xq) with law Pq (defined in (jlip ) 
independent of B . Then, the nonlinear stochastic differential equation 

i) Xt = XQ + V^uJq l{s>^}dBs + jQK{p){s,Xs)l{s>T}ds, 

^t = h + JqVK{p){s, Xs)<^>sl{s>r}ds, and (28) 



It 



Hi) the law P o/(r, X, <I>) belongs to V^s and Pt{dx) = p{t,x)dx 



has at most one pathwise solution. Moreover, its law is a solution of (MP). 
c) If the condition 



^ (||w^ol|p + r|||g|||o,p,0) <ro(p) 

is satisfied, where Tq{p) > is the constant provided by Theorem \3.5l then a solution 
to (MP) in VJ^s exists, and there is strong existence for ([28]) in [0,r]. 

"'2 

The proof of Theorem 14. 21 requires some prehminary facts about a scalar problem implicitly 
included in the vectorial problem (MP). 

4.1 A nonlinear Fokker-Planck equation with external field associated 
with the 3d-vortex equation 

For any Q € , we denote by Qt the sub-probability measure on defined for scalar 
functions by 

Qtif) = i?«(/(Xi)l|,<i}), (29) 

where (r, X) are the two first marginal of the canonical process (r, X, $) in Ct- Obviously, 
for Q S we have 

Qt-^Qt-^ Qt, 

and we shall denote 

kf{x) := ^(x), (30) 
dQt 

and, indeed, 

,Q, , EQ{'^th{T,Xo)l^,^t}\Xt = x) ^ 

h (^) = Q{r<t\Xt~=x) l{Q(-<*l^*-)>0}- 

If Qf{dx) has a density p'^{t,x) with respect to Lebesgue measure, we shall respectively 
denote by fP{t,x) and p'^{t,x) the densities of Qt and Qt- Then, clearly one has 

pit,x) =kfix)p^it,x). 

Remark 4.3 If Q € VbiCx) is such that Qt is absolutely continuous for all t G [0,T], the 
existence of a joint measurable version of {t,x) ^ p^{t,x) is standard by continuity of Xt 
under Q°. We always work with such a version. Moreover, there exists measurable versions 
of {t,x) 1-^ p^{t,x) and {t,x) i— > p^{t,x). This can be seen by Lebesgue derivation (see e.g. 
Theorem 3.22 in f^), taking 5 — > in the quotients 

Q{T<t,XteB{x,6)) EQ{'^tHT,Xo)l{,^t},Xt e B{x,6)) 

QiXtGB{x,6)) QiXteBix,S)) 
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and using the previous relation between p^(t,x) and (here, B{x,6) is the open ball of 
radius r centered at x). 

Lemma 4.4 Assume that (MP) has a solution P € VbiCx) such that has a density for 
each t £ [0, T] and such that condition (I15|) holds. Then: 

i) the couple {p,p) satisfies the weak evolution equation 



fit,y)p{i^y)dy= / fiO,y)wo{y)dy + 



fis,y)g{s,y)dyds 



+ 







ds 



s, y) + u/\f{s, y) + K(p)(s, y)V/(s, y) 



p{s,y) dy ds, (31) 



for all f G C^'^, where wq and g were defined in (llOp . 
a) p is moreover a solution of the mild equation in [0, T], 

p{t,x) =G'^*wo{x) + [ G'^^,*g{s,-){x)d. 
Jo 

dGt 



rt 3 



(32) 



{x -y)K{kp)j{s,y)p{s,y) dy ds, 



where the multiple integral converges absolutely and k is the function k^ defined in 



Proof, i) From the definition of (MP), and the fact that l{t>r} is J'o measurable, we see 
that the expectation of the expression 



f{t, Xt)l{t>T} - f{T , -'^o)l{t>r} 



-^{s, Xs) + uAfis, Xs)ds + K{p){s, Xs)Vf{s, X,) 



l{s>T}dsl{t>T} 



vanishes. On the other hand, condition (fT5]) obviously implies that 

/ \K.{p){t, x)\p{t, x)dxdt < oo, 

J[0,T]xR3 

which allows us to apply Fubini's theorem after taking expectation and conclude the asser- 
tion. 

iijFix^ GV and t € [0, T] and take in ([311) the C^^'^-function ft : [0,t] x ^ ^^^^^ 
by ft{s, y) = Gt_g * il^{y) (which solves the backward heat equation on [0, t\ x with final 
condition f{t,y) = tp{y)). By Lemma |3. II and condition (jlSp . it is not hard to check that 

^{p)j(.s,y)\\tlj{x)\p{s,y)dx dy ds < oo. 

By Fubini's theorem we easily conclude. [] 



r r y dG^ 
Jo im3)2f-f dyj ^ ^ 
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Consider now a fixed but arbitrary function k : [0, T] x ^ M'^ of class L°° {[0,T], Lf), 
and formally define an operator on functions rj,v ^ A4eas^ by 

h^{r],u){t,x) = V/ -—^{x-y)K{ku)j{s,y)r]{s,y)dyds. 
Jo jr^^M3 oyj 

If now Fq^p^t, Fi,p,T, -^o,r,{T;p) i*'i,r,(T;p) denote the scalar-function analogues of the 
spaces F defined in Section 3, we have 

Remark 4.5 For each p £ [1, oo] (resp. each p G [1, oo] and r > p), the mapping rj ^ krj 
is continuous from Fq^p^t to Fo,p,T (resp. from -Fo,r,(T;p) to FQ^r,{T;p))- 

Writing now 

7o(t, x) := * wo{x) + / Gt^^ * g(s, ■){x)ds, 

Jo 

we can state the following properties of the scalar equation (j32p . 

Proposition 4.6 Assume (Hi) and (Hp) to^/z p G (|,3), and let k G L°°([0, T], Lf") be a 
fixed but arbitrary function. 

i) For each r G [p, oo), we have 

70 G ^0,r,(T;p), "^^^/l lllfo lllo,r,(T;p) < C{r,p)\\wQ\\p + T|||g|||o,p,r 

/or some finite constant C{r,p) > 0. 

ii) Suppose that |<p<3, p < I < rnin{-^^,3} and ^ < I' < gz^- T/ien, i/iere exists 
a finite constant Cq{1, l';p) not depending on T > such that for all ?/, G -^o,i,(T;p)) 

\l^Hr],'y)\lo,l',{T;p) < Co(/,^';p)r^~^|||r/|||o,/,(T;p) III l^lllo,/,{T;p) • 

Hi) The mild Fokker-Planck equation with external field i32\) has at most one solution 
p G Fq^p t for each T > 0. 

iv) If p £ -^o,p,T is a solution of /f^). then p G -Fo,r,(r;p) for all r G [p, oo) with |||p|||o,r,(r;p) ^ 
C{T,p,r, |||p|||o,p,t) < oo. 

v) We deduce that for all I G [3^,oo), K.(kp) G F-, , 3p y 

Proof. Part i ) follows from Lemma 13.11 in a similar way as part of Lemma 13.41 We 
notice that the restriction on r in the latter was needed only to ensure that the derivative 
of time integral was convergent, and so it is not needed here. Thanks to Remark 14. 5^ part 
ii) is similar to part ii) of Proposition 3.1 in [TO] . 

Prom the previous parts, equation (j32p admits the abstract formulation in Fq^p^t 

/5 = 7o + b'=(p,p). 

Then, the arguments yielding parts i) oi Theorems 13.51 and 13.71 also provide the assertions 
of parts Hi) and iv)., repectively. For part uj, we notice that from iv), kp G Fo^j.^(r;p) holds 
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for all r G \p, oo[. Thus, if we take l> q := and set r := ( -j- + i) ^ then one has r > p, 
and so Lemma 13.31 i ) implies that 

sup t2Vp r^||K(A;p)(t, •)||; = sup ; 
te[o,T] te[o,T] 



sup t2^p r)||K(A;p)(t, Oil; = sup ta^? '^||K(A:/5)(t, ^Ib < oo. 



This shows that K(A;/5) G 'Po,i,{T;q)- We conclude that K(A;/5) G F;^ ; (-p.^), noting that since 

kp G Fo^/^(T;p) for all / > q, Lemma [HTU] i) implies that ^^^^^^ G ^o,l,{T;p) for all A: = 1,2,3. 
In other words, 



sup 

t£[0,T] 



dK{kp){t,' 



dxi 



sup t2^2S l) 



I te[o,T] 



dK{kp){t,' 



dxi 



< oo, 



□ 



which is the required estimate. 

4.2 Uniqueness in law and pathwise uniqueness 

We need the following version of Gronwall's lemma: 

Lemma 4.7 Let g and k be positive functions on [0, T], such that k{s)ds < oo, g is 
bounded, and 

g{t)<C+ / g{s)k{s)ds for all t e [0,T]. 
Jo 

Then, we have 



g{t) < Cexp f k{s)ds for all t G [0,r]. 
Jo 



We are ready to prove parts a) and b) in Theorem 14.21 

Proof. Let P G VTa be a solution of (MP). Since p G -Fo,i,T H i^o.p.Ti by interpolation we 



have p G 3rp. By Lemma [313] i) we deduce that (fT5]) holds. Moreover, by Lemma ii). 
Proposition 14.61 iv) and Lemma [STBl i). we have that VK(p) G -Fo,3, (T;p) smd, consequently, 
condition (jl6p also holds. By Lemma 12.61 we deduce that p is a weak solution of the vortex 
equation and, since /cf is bounded, we have p G Fo,p,T- 

We now need to prove that the latter implies that p G Fo,p,T is uniquely determined. By 
Theorem 13.51 a) this will follow by checking that p is also mild solution. For fixed "0 £ ^3 
and t G [0,T], define ft ■ [0,t] x ^ by ft(s,y) = G'^_, * V'(y), which is a function 
of class (C^'^)3 that solves the backward heat equation on [0, t] x M.^ with final condition 
i{t,y) = ip{y). One can thus take ^t in the weak vortex equation and, thanks to conditions 
(I15p and ()16|) . apply Fubini's theorem to deduce (since -0 £ ^3 is arbitrary) that 



p{t,x) 



Jo ~[ JR" 



{x - y)[K(p)j (s, y)p(s, y)] 



+ GU{x-y)[p,{s,y)^-^{s,y)] 



dy ds. 
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Since p is divergence-free, to see that p solves the mild equation it is enough to justify an 
integration by parts of the last term in the previous equation. We cannot do that at this 
point since we cannot ensure enough (Sobolev) regularity of p. But noting that for q = 
one has 1 < < |, we see that the function p = k^p belongs to Fo,g*,T by interpolation. 
On the other hand, one has — •)K(p)(s, •) G W^''^ thanks to Proposition l4.6l f ). Since 

by hypothesis, div p{s) = in the distribution sense, the fact that p{s) S and a density 
argument allow us to check that 

E / ~Pj{s,y)^[GU{x-y)K{p){s,y)]dy = 

for all s g]0, T]. Thus, w := /5 is the unique solution of p8]) in Fo,p,T- 
Now, by a standard argument using the semi-martingale decomposition of the coordinate 
processes X* and their products X^X^, we obtain that the martingale part of f{t,Xt) 
in (MP) is given by the stochastic integral f^V f{s, Xs)l{s>T}dBs, with respect to a 
Brownian motion B defined on some extension of the canonical space. From this and the 
previously established uniqueness of p, P is the law of a weak solution of the stochastic 
differential equation 

i) Xt = Xo + V2^ /o l{s>r}dBs + /o K(w)(s, {s>T}ds, 

(33) 

ii) <^t = h + f^VK{w){s,Xs)^sl{s>r}ds. 

Since the equation (j33p is linear in the sense of McKean, to conclude uniqueness in law it 
is enough to prove pathwise uniqueness for it. This is done first for X and then for both 
with help of the estimate on ||VK(w)(t)||j^ in Theorem 13.71 and Gronwall's lemma. [] 



4.3 Pathwise convergence of the molhfied processes and strong existence 
for small time 

To prove part cj of Theorem 14. 2^ we shall construct a solution via approximation by 
nonlinear martingale problems with regular drift terms K^(w^) and VK^(w^), where 
G -^o,p,Tn-Fo,i,T is given by Theorem l3.5[ Our arguments will improve the ones developed 
|10j . and yield a pathwise approximation result and an explicit rate. 
If q = -0^, Holder's inequality and the properties of K imply that that for all t £ [0,T], 

||K=(w^)(t,-)||oo < C||(^,||,*|||K(w^)|||o,,,T < Cye\\g4^'\h,p,T- 

Similarly, one has ||VK^(w^(t))||oo < C||V99e||g* |||w^|||o.p.T and analogous estimates hold for 
all derivatives. Thus for each e > 0, the function (s,y) •— > K^(w^)(s,y) is bounded and 
continuous in y G M^, and has infinitely many derivatives in y G M^, which are uniformly 
bounded in [0,r] x M^. 

We fix now the time interval [0, T] given by Theorem 14.21 It will be useful to consider in 
this section the stochastic flow 

es,t{x) = x + V2^{Bt - Bs) + fK%w'){6,C,e{x))de, for all t G [s,T], (34) 

J s 

which has a version that is continuously differentiable in x for all {s, t) thanks to the 
regularity properties of K^(w^) (cf. Kunita [17J). 
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We also consider the strong solution of the stochastic differential equation in [0, T] 

Xf = Xo + V2^[ l{s>r}dBs+ I K^(w^)(s,Xf)l|,>,|ds 
^ Jo Jo 

Jo 



(35) 



where (r, Xq) is independent of B. We denote by the joint law of (r, X^ , $^) and observe 
that G V]^ . Since Xf = Xq for all t < r, we have that 

= ir,t{Xo)'^{t>T} +-'^0l{t<r}- 



Denoting by G^(s, x; t, y), (s, x, t, y) € 



p2\2 



,s<t the density of ^Iti^) (which is a 



continuous function of {s,x,t,y), see [13j), and conditioning with respect to {t,Xq), we 
obtain for bounded and measurable functions / that 



f{y)G^s,x;y,t)dyPo{ds,dx) + 



^(IR3)2 

f{x)'Wo{x)dx + 
+ 



f{x)Poids,dx), 



JR3 



f{y)G''{s,x;t,y)dy 



g(s, x)dx ds 



T 



/(x)g(s, x)dxds. 



Consequently, Xf has a (bi- measurable) family of densities that we denote by . Observe 
that one has p'^{t) G for all t E [0,T] from the assumption on wq and g and standard 
Gaussian bounds for G^(s, x; t, y). 

The functions and correspond to the densities of, respectively, the sub-probability 
measure and the vectorial measure 

f^E [f{i^^t{Xr))l{t>r}\ and f ^ ii; [f(e,i(X,))V,e^,,(XO/i(r,Xo)l|t>,}] . 

They are bi- measurable by similar arguments as in Remark 14.31 and we have p^ (t) G 
and p%t) e Lf . 

The assumptions on (p ensure the following estimate concerning the approximations of 
the Dirac mass, see Lemma 4.4 in Raviart 



Lemma 4.8 Let if be a cutoff function of or order 1. Then, for all v G W^'^ and r G [1, oo], 
one has 



\v — ipe * v\\r < Ce 



i=l 



dv 



dxi 



We deduce the following result: 

Lemma 4.9 i) We have pF = w'' and, consequently, 

sup |||p^|||o,p,r < oo and sup |||/)^|||o,p,T < oo. 

e>0 ' ' £>0 



(36) 
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a) If (f is a cutoff function of order 1, the we have that 



sup 2 _ w(t)||p < C{T)e 

iG[0,T] 



for some finite constant C{T). 



Proof, i) Since E [J^ \K^{w^){t, Xf)\dt\ and E [J^ |VK"(w^)(t, Xf j are finite, we 



can follow the lines of Lemma 2.6 and use Remark 2.4 to see that for all f G (C^' )3([0,r] 



i{0,y)'Wo{y)dy + 



di 

— (s,y) +i/Af(s,y) 







Jo Jr3 

+Vf (s, y)K^(w^)(s, y) + f (s, y)VK^(w^)(s, y) 



Hs,y)^{s,y)dy ds 



P''{s,y) dy ds. 



(37) 



On the other hand, the regularity properties of the stochastic flow ()34p imply that for all 
<j) £ V and 9 g]0, T], the Cauchy problem 

f/(s,y) + z.A/(s,y) + K-(w-)(s,y)V/(s,y) =0, {s,y) e [0,e[xR^ 

(38) 

/(0,y) =,^(y). 

has a unique solution / that belongs to Cl'^{[0,9] x M^) (see Lemma 4.3 in [10]). One 
can thus use the function f = V/ in equation (f37|) . and after simple computations obtain, 
thanks to the null divergence of wq and g(s, •), that 



V</.(y)p(")(t,y)dy 



V 



ds 



(s, y) + ivA/(s, y) + K^(w^)(s, y)V/(s, y) 



p(")(s,y) ds = 



for all (j) £ v. Thus, div p^{t) = 0, and we can adapt the arguments of Section 4.2 to 
conclude that solves the linear mild equation 



V = Wo + B^(v, 



V e Fo,pT- 



(39) 



Since uniqueness for (j39p holds (by similar arguments as for the nonlinear version), and 
also solves the equation, we conclude that = w^. The asserted uniform bound for is 
thus granted by Theorem 13. 5[ To obtain the uniform bound for p^, we take L*' norm to 
(139p . and follow the arguments of the proof of Theorem 13.51 i), to get that 



p < |||wo|||o,p,T + C'|||w^|||o,p,T / (t-s) \\p%s)\\p ds, 



from where conclude applying Gronwall's lemma similarly. 

ii) By an iterative argument as in the proof of Theorem 13.51 we get that 

\\p'{t)-^{t)\\p <C t a{t-s) ||K^(w)(s)-K(w)(s)||g 
Jo 



+ C{T) f 
Jo 



(40) 



|p^(s) - w(s)||, ds, 
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where a{s) = Xlfc^i^ ^) = 1 ~ j" and N{p) = [9q ^\ + 1. Integrating in time and 
using Gronwall's lemma, Theorem 13.71 ij and Lemma |4.8^ we obtain that for all 9 G [0,T], 



/ / a{t- s) \\K'{w){s) -K{w){s)\\q ds dt 
Jo Jo 



\p'{t) -w{t)\\pdt <C 



<Ce 



t^^^--v^'2dt = eC{T). 

k=l 



Injecting the latter in (j40p . we obtain 

\\p%t)-^{t)\\p <eC{T) + C f a{t-s) ||K^(w)(s) - K(w)(s)||, ds 

Jo 



<eC{T) + Ct^'^e 



and the conclusion follows. 



□ 



The Proof of Theorem 14.21 c) will be completed with the following result, which moreover 
establishes the strong pathwise convergence of the nonlinear processes (X^,<I>^) as e ^ 0. 
We are inspired here in ideas introduced in [12], but we need a finer use of analytical 
properties, as we shall improve the rate of , 6 € (0, 1), that was obtained therein for a 
particular choice of kernel. Further difficulties also will arise because of the additional (and 
more singular) drift term of the "vortex stretching processes", proper to dimension 3. 

Proposition 4.10 Let if be a cutoff of order 1 and defined in terms of ip as be- 
fore. Then, as e goes to 0, the family of processes {X^ — Xo,^^),£ > is Cauchy in 
the Banach space of continuous processes (Y, ^) with values in x R'^®'^ with finite norm 
-E'(sup4gjo,T] \^t \ + \^t\)- Moreover, one has 



E sup \Xt - Xt\ + - $f I < C{T)e, 



where {X,^) is a solution of the nonlinear s.d.e. 

Proof. We observe that the substraction of Xq is only needed to avoid a moment-type 
assumption on Xq. Let e > e' > 0. We have 



ELuv\Xl-Xt\] < f E (K^(w^)(s,Xf)-K^'(w^)(s,Xf))l|,>,} 

\s<t J Jo 

+ f E (K-'(w^)(.,Xf)-K-'(w^')(.,Xf))l|,>,} 
Jo 

+ f E (K^'(w^')(.,XD-K^'(w^')(.,xf))l|,>,| 
Jo 



ds 
ds 
ds 



(41) 



The third term on the r.h.s. of (|4ip is bounded thanks to Theorem 13.71 Hi) by 



c[ s'^'^^p'-'^E (snp\Xl - Xl'\] ds, 
Jo \e<s J 
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for any fixed r G (3, -0^)- Writing q = and q* for its Holder conjugate, and using 



Lemma 13.31 and Lemma 14.91 ii), we bound the second term by 

[ ||K^'(w^)(.) - K^'i^^')is)Up%s)yds < C{T)e. 
Jo 

We have used the fact that sup^^g lll/5'^lllo,g*,T < cjo by interpolation since g* < | < p. By 
similar arguments, the first term on the r.h.s. of ()4ip can be bounded above by 



Jo 



|K-'(w=)(s) - K^^n{s)\UW{s)\U*ds < C{T)e. 



Bringing all together and using Gronwall's lemma we deduce that 



E sup \Xl -Xt\\< C{T)e. 



(42) 



Now, notice that Gronwall's lemma and Theorem 13.71 in) imply that the processes are 
bounded uniformly in t £ [0,T],e > and randomness. Therefore, we have 



(43) 



ii;fsup|cl>^-c^f|') <C ['e (VK-(w-)(5,Xf)-VK-'(w^)(s,Xf))l|,>,| 

\s<t J Jo 

+C C E (VK-'(w^)(s,Xf)-VK^'(w-')(s,Xf))l|,>,| 
Jo 

+C ( E (VK-'(w-')(s,Xf)-VK^'(w^')(s,xf))l|,>,} 
Jo 

+C [ El VK^'(w"')(s,xf ) sup|^>^ - <^l'\] ds 
Jo \ e<s J 

By Theorem 13.71 Hi), for fixed r S (3, q) the last term in the r.h.s. of (j43p is bounded by 

ft 



ds 
ds 
ds 



C I I sup - $|'| I ds, 

'0 \e<s I 



and the third one is by 



1 'i(l-l) 



/o 



ds < C{T)e, 



using also the previous estimates on E\Xl — XI \ . The first term in (j43|) is upper bounded 
by 



C [ W{s)\\p' ||VK=(w^)(s) - VK''{w'){s)\\pds. 
Jo 



(44) 



If p > 2, then we have p* < 2 and so by (j36p and interpolation, we deduce that (j44|) is 
bounded by 

C\W\lo,r,T r ||VK((^, * w^)(s) - VK(w^)||p + llVK(w^) - VK((^,, * ^'){s)\\pds 
Jo 

< CTe. 
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This last inequality is obtained by Lemma [3.6l i). LemmaHSl Lemma [4.91 i) and the uniform 
boundedness of (w'^)£>o in Fi^p^T- If now | < p < 2, then we have 3 > p* > 2 > p and by 
similar steps as in the previous case p > 2, we can upper bound (pil) by 

C|||/5^lllo,p*,(T;p) r s-|^-p-^^\\VK^^^^'){s) - VK^'(w^)(s)||pd. 
Jo 

< esup|/5 ||o,p*,(T;p) / s 2> p*'s '^ds 

<5>0 JO 

< eC{T) sup |||/|||o,p*,(T;p)- 

<5>0 

We have used here Lemma HTHl the fact that (w'^)e>o is uniformly bounded in ^i^p^r and 
that ~|(| — pr) — ^ > —1 since p > |. The fact that the supremum in the previous 
estimate is finite, is seen in the same way as part vi) of Proposition 14.61 namely by an 
iterative argument using the mild equation (similar as therein) satisfied by p^, starting 
from the uniform bound in Lemma 14.91 i). 

Thus, we have shown that the first term in the r.h.s. of (j43p is bounded by a constant 
times e. Let us now tackle the second term in the r.h.s. of (j43|) . This is bounded by 

C r \\p'{s)\\p* II VK^' (w^) (s) - VK^' (w^') (s) llpds 
Jo 

<C [ \\p'{s)\\p*\\w'{s)-w''{s)\\pds, (45) 
Jo 



thanks to Lemma 13.61 By Lemma 14.91 ii ) we can upper bound (j45p respectively by 

/o 



Ce I s^-^ds = eC{T), 



in the case p > 2, or by 

rT 



f 3 / 1 1 N 1 3 

Jo 



in the case p < 2, where the constants are finite since p > |- 
Consequently, we have an estimate of the form 

E fsup 1$^ -<^l'\) <Ce + C [ s~^-^^p-r'^E (sup |$^ -^'\] ds, 
\s<t J Jo \e<s J 

for each fixed r G (3, g), and Gronwall's lemma yields 



E { sup|$^ < C{T)e (46) 



for all e > e' > 0. 

Estimates (j42p and (j46p thus show that (X^ — Xq, <1>^) is a Cauchy sequence in the Banach 
space of continuous processes {Y, ^') with values in R'^xR'^'^^ and finite norm i?(supjg[o,T] 1^1 + 
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I^'tl). Write the limit in the form {X — Xq, $), for a continuous process {X, $) and define 
£} and by the relations: 

Xt = Xo + V2^ [ l{s>r}dBs+ ! K(w)(s,X,)l|,>,|ds + f:/, 

l-t ''^ (47) 
^t = h+ / VK(w)(s,X,)$,l|,>,}ds + f2_ 
Jo 

Comparing (X, $) and (X^,$^), and using similar estimates as so far in this proof, but 
with instead of e' (and w instead of ), we get that (X, $) satisfies (|47p with £1 = 0, 
i = 1,2. Since that is a linear s.d.e (in McKean's sense), the proof that (X, $) is the 
asserted nonlinear process will be achieved by checking that for all bounded Lipschitz 
function f : ^ R^, one has 

^(f(Xt)$t/i(r,Xo)l|,>,|) = / f(x)w(t,x)dx. 

This follows from the facts that 

E{i{Xl)^h{T,Xo)l{,>r}) = I f(x)w^(t,x)(ix 

and 

|E(f(Xi)$i/i(T,Xo)l|,>,}) - i?(f(Xf)<Df/i(r,Xo)l|,>,})| 

< (||<i>||Loo([o,T]xn) + l)\\h\\MLipE{\Xt - Xl\ + - ^\) 

<C\nLipe. (48) 
□ 

Remark 4.11 a) By Lemma i), the process (X^,<I>^) defined in (135p is a solution 
in [0, r] of the nonlinear s.d.e. 

i) Xf = Xo + l{s>r}dBs + /o*K^(/5^)(s,Xf)l|,>,|ds, 

ii) $f = /3 + /oVK^(p^)(s,Xf)$^l|,>,|ds, and (49) 



ni) t/ie law o/ (r, X^, <I>^) belongs to Vj^^ and P^{dx) = ff{t,x)dx. 



b) It is also possible to associate a unique pathwise solution of (128p with any solution 
w € Fo,p,T n Fo,i,T of the mild vortex equation (i.e. not necessarily the one given by 
Theorem \3. 5\) . This can be done by an approximation argument similar to the previous 
one, but considering linear processes in the sense of McKean (with drift terms K^(w) 
and VK^(w)J instead of the processes (j35]) . 



c) Denote now by Wt the Wasserstein distance in V{Ct) associated with the metric in 
Ct := [0,r] X C([0,r],M3 X 



d{{0,y,^),{ri,x,(t))) =\9 -r]\ + sup (mm{|x(t) - y{t)\,l} 

te[o,T] 

+min{\i;{t)-^{t)\,l]). 
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Then, the previous proof states that 

WTiP',P)<C{T)e, 

where P is the law of the nonlinear process (j28l) . 

d) By the regularity results of Section 3, one can prove in a similar way as in Corollary 
4-3 of llpl that the stochastic flow ([6]) is of class , in spite of the fact that u and 
Vu are singular at t = 0. Thus, the identity ([7]) holds. 



5 The stochastic vortex method 

We consider first a cutoffed and mollified version of the vortex equation with external field 
(l2|) which extends the McKean-Vlasov model studied in [lOj when g = 0. 
Denote by Mg the sup- norm of i^^ on and by a Lipschitz constant for it, which 
respectively behave like p- and ^ when e << 1. Notice that div = [div K) = 0. 
For it! > 0, we denote by XR '■ M^®^ M'^®^ a Lipschitz continuous truncation function 
such that |X-r(0)I ^ P- We may and shall assume that XR has Lipschitz constant less than 
or equal to 1. 

Consider now a filtered probability space endowed with an adapted standard 3-dimensional 
Brownian motion B and with a [0,T] x M^-valued random variable {t,Xq) independent of 
B and with law Pq. 

Theorem 5.1 There is existence and uniqueness (pathwise and in law) for the nonlinear 
process with random space-time births, nonlinear in the sense of McKean: 



^ Jo Jo 

a>^'^ = /3+ / Vu''%,Xl''')xRm'')l{s>r}ds 

Jo 



ds 



(50) 



with 

u''%, x)=E [K,{x - '^O A XR{^'s'')h{T, Xo)l|,>,|] . (51) 

The proof is based in the classic contraction argument of Sznitmann [25] and is not hard 
to obtain by combining elements of Theorems 5.1 in [lOJ and Theorem 3.1 in |12j . 

Consider next a probability space endowed with a sequence (i?*)jgN of independent 3- 
dimensional Brownian motions, and a sequence of independent random variables (r*, XQ)jgp^ 
with law Pq and independent of the Brownian motions. For each n G N and R,£ > we 
define the following system of interacting particles: 

X.,e,R,n + ^ /■*! di?^ 

Jo 

+ f -Y.K,{Xl'^^^'- - X^'^'") A Xi^($r'^•")/^(r^X^)l^,>,.,,,}ds 

^i e,R,n (52) 



+ 



XR{^t''nHs>r^,r^}ds, 
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for i = 1 . . . n, and with VK{y) Az = Vy{K{y) A z) for y, z & M'^, y 7^ 0. Pathwise existence 
and uniqueness can be proved by adapting standard arguments, thanks to the Lipschitz 
continuity of the coefficients. 

In the same probabihty space, we also consider the sequence 

XY-'' = Xl + V2^ j l|,>,.}dS:+ fu''%,Xl''''')l{,y,.yds 

t •'^ -^0 ,iGN (53) 

= /3+ / W^R{s,Xl'^^^)xR{^T'')Hs>r^}ds 

Jo 

of independent copies of (j50p . Their common law in Ct is denoted by P^'^, and we write 
h := ||wo||i + ||g||i,r- Recall that XR is Lipschitz-continuous function, bounded by i? > 
and with Lipschitz constant less than or equal to 1. It is not hard to adapt the proof of 
Theorem 5.2 in [lOj to get the following: 

Theorem 5.2 For e > sufficiently small and all R > 0, we have 



E 

for all i < n, where 



{\ ^i.e.R.n ^i,e,R\ , \^i,e,R,n ii,£,R\\ 



<^C{e,R,h,T) (54) 
'n 



C(e, R, h, T) = Cie(l + RhT){RhT) eicp{C2e-^hT{R + l){h + RT)} 
for some positive constants Ci,C2 independent of R, £,T and h. 

Let us now make the assumptions of Theorem 13. 5( and consider, in the corresponding time 
interval [0, T], independent copies (X*'^, <I>*'^) and (X,<I>*) of the processes (f28|) and (fi9|) 
constructed on the given data {Xq,t^ , B^), i G N. 

Recall again that the uniform bound of Theorem 13.71 Hi ) and Gronwall's lemma imply that 
the processes are uniformly bounded, say 

sup |$f(w)| < i?o(T,wo), 

t£[o,T],£>o,uj<=n 

for some finite positive constant Ro{T, wq). Thus, for any R> Ro, one has for all t £ [0, T] 
that 

{xi'^<^r) = {xr,XRK'')). 

Consequently, (X*'^, $*'^) is a pathwise solution in [0, T] of equation ([55]) . and so we conclude 
that 

almost surely. Bringing all together, we obtain the following pathwise approximation result: 

Theorem 5.3 Assume that (Hi) and (Hp) hold with p G (|,3) and that the hypothesis of 
Theorem \3.5\ i) is satisfied. Let Kg, he defined as in (|19p . with ip a cutoff function of order 
1 and write h = ||wo||i + ||g||i,T- Let furthermore R > Ro{T,wq) and 

£n = (cQ,lnn)"9, 
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with 



0<Ca< a{C2hT{R + l){h + RT))-\ 



for some alpha a G (0, i). Then, we have for all i < n, 



E 



sup -Xi\ + 1$^'^"'^'" - 

telo,T] ^ 



< C{T,wo,g,a) 



+ 1 

n2 "(Inn) 9 (lnn)9 



, (55) 



where (X,^) is the unique pathwise solution of (|28p . and the constant C{T,wo,g,a) de- 
pends on the data wq and g only through the quantities \\wq\\p, |||g|||o,p,T cind \\wo\\i + ||g||i,T- 

Remark 5.4 i) The rate at which the second term in the right hand side of ()55p goes 
to is exactly that of e = En- The logarithmic order of latter was needed to make the 
upper bound in Theorem \5.2\ go to with n, which then happens at an algebraic rate. 
The global rate is therefore conditioned by the techniques used in the proof of Theorem 
15-4 see ilO^ for details. Under additional regularity assumptions, it is possible by 
analytic arguments to slightly improve the convergence rate (see the discussion at the 
end). An attempt for a more substantial improvement should however exploit specific 
features of the interaction at the level of the particle systems. 

a) The previous result implies as usual that Wt (/ai(;(X*'^''^'", <I>*'^'^'"'), P) goes to at 
least that fast, and that (with the obviously extended meaning ofWr) 

Wt (law (^(X^'^'^'", $i'^.«'"), . . . , (x^'^'^'", ^k,s,R,n^^'^ ^ p®fc^ < ^^^^ 

where 5n stands for the quantity in the right hand side of (j55p . 

We deduce the convergence at the level of empirical processes: 

Corollary 5.5 Under the assumptions of Theorem \5.3\ the family 
(/i"'^"'^)o<t<T of M.^ -weighted empirical measures on M.^ 



^n,e^,R ,= l.^S^.,.^,n,r. ■ (xij($r"'^'")/io(r,X^)) 1|,>,} 
i=l 



converges in probability to {w{t,x)dx)Q<t<T in the space 

C([0, r], 7W3(]R^)), where A^3(M^) denotes the space of finite M.^-valued measures on 
endowed with the weak topology. Moreover, we have 



sup E 

te[o,T],||f||i„p<i 



(^r"'''-w(t),f) 



< C 



1 



+ 



1 



- + 



1 



n2 ""(Inn) 9 (Inn)! 

where ||f|| Lip is the usual norm in the space of bounded Lipshitz continuous functions f : 
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Proof. It is enough to prove the bound for Lipshitz bounded functions. For such a function 



~n,£n,R £\ 



(w(t),f) 



< 



+ 



~n,en,R 



1 " 



1 



(56) 



n 



i=l 



Ct 



f (y(i)) A XRm))H9, xmP^'-'^id9, dy, 
+ Kw-(t)-w(t),f)| 



with P^"'^ = P^" = law{T, X*'^"'-^, (I)*'^"''^^). The independence of the processes (r*, X*'^'"^, $*'="'^) 
i G N and the definition of h imply that the expectation of the second term in the right 
hand side of ([56]) is bounded by ■^2\\i\\LipRh, where h = (||tfo||i + ||g||i,T)- We use the 
latter and estimate in Theorem 15.21 to bound the first term, and get that 



E 



(/2,"'""'^ - w(t), f) < \\nLip{R + l)/i^C(e„,, R, h, T) 



1 



n 



+ 



2||f lliipii/l 



'n 



+ |(w^"-w(t),f) 



The last term being equal to the first term in (08]), the conclusion follows. 



□ 



Remark 5.6 In the case g = 0, Philipowski ]22^ obtained a similar approximation result 
for the vorticity field (for a simpler particle system) under the additional assumption that 
the test function f belongs to . 

Finally, we establish an approximation result with convergence rate for the velocity field. 
To that end we will need to strengthen to convergence of to w already shown. We shall 
need the following 

Lemma 5.7 For each p G (|,p), there is a constant C{T,p) such that 



sup 



t^\\Vw%t) - Vw(t)||p < C(T,p)e. 



Proof. We need p G (|;3) in order to dispose from a integrable (in time) bound for 
||Z'^K''(w^)(t)||_3£_, which we do not have for p = p. Indeed, for any p in that interval 



3-p 

3p 



we have q := G (3, g^), and so by Theorem 13.71 i) and Lemma 13.61 we have for 
k,j, i = 1, 2, 3 that 



with 



-(---) 
sup P 1 

te[o,r],£>o 
,1 _ 3^1 _ 1) = 



dxi 



+ sup t2^2yp g> 

q te[0,T],e>0 



dxjdxi 



< oo, 



-1 + |(i — |) > —1- Let us now check that one has 

sup |||w^||i^p^T' < OO. 
e>0 



(57) 



(58) 
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This is not immediate, since T > given by Theorem 13.51 was determined by the norm of 
wq and of the operator in the spaces corresponding to the parameter p > p. We wih 
prove ()58p using continuity properties of the operators B^. It follows from Proposition 3.1 
in) in 



|10j that for | < r < 3 and < r' < r, one has 



sup |||B^(v, v)|||iy_T < Cr, 

e>0 



I ^2 



(59) 



for some finite constant Cr^r'iT)- From this, we deduce that G Fi^p^T, with a uniform 
(in e) bound, by the following iterative procedure. Define a real sequence by ro = p, 
r„+i = , and notice that it is increasingly convergent to 3. We can thus take G N 
such that r]\i < p < r^v+i. The function s i-^ being increasing on [0,6], we then have 



3p ^ 3rjv^i 
6-p — 6-rjv+i 



r]\f. By ([59]) with r = p and r' = r^, we see that B'^(w'^, w^) G ¥i^ri^,T-, and 
since also wq G Fi^^jv.t holds by Lemma [3^ i) (taking in the place of p and r therein), 
we get that G Fi^rjv,T, with a bound in that space that is uniform in e. We repeat the 
previous arguments with r = and r' = = ^^Af-i and get that w"^ G Fi^rj^_i,T, with 

a bound that is a uniform in e. Continuing — 1 times this scheme we get ()58p . 
We now take derivatives in the mild vortex equation with e > (as justified in the proof 
of Proposition 3.1 in [TO]), 



dxi 



{t,x) 



QQu rt 

— ^(x - y){wo)k{y)dy + 

OXi 



dxi 



-(x - y) 







5G| 
dxi 



{x - y)g{0,y)dy ds 



K^(w^),(s,y) 
- w^(s,y) 



^w|(s, j/) 

9yj 

dK'{^%{s,y) 



dy ds, 



for k = 1,2, 3. Notice now that, thanks to the estimates (j58p . Lemma [4. 91 ii) also holds with 
p replaced by p. By estimates as those in the proof of Theorem 13.51 i ) and using Lemma 
14.91 ii ) and estimates ()57p and (|58p , we then that have 

||Vw^(t) - Vw(t)||p < C [\t - s)-*s-5 [||w^(s) - w{s)\\p 

Jo 

+ \\K%v,){s)-K{^){s)\U]ds 

+ C [ (t-s)-i[||Vw^(s)+Vw(s)||p 
Jo 

+ ||VK"(w)(s) - VK(w(s)||,-]ds 

3 /* ^ 3 3 3 

<Cet^~'p+Cel {t - s)~'^s~^^^~^ds 
Jo 

+ C [ (t-s)"^||Vw^(s) - Vw(s)||pds 
Jo 

<Cet~^+C [ (t-s)~^||Vw"(s) - Vw(s)||pds. 
Jo 

Iterating the latter sufficiently many times (using the identity quoted in the proof of The- 
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orem [33]) i), we obtain that 

||Vw^(t)- Vw(t)||p<Ce(ri + l) + C7(r) /" ||Vw"(s) - Vw(s)||pds. (60) 

JO 

Integrating (j60p in time and using Gronwall's lemma, and then inserting the obtained bound 
in the r.h.s. of ()60p . we obtain 



||Vw^(t) - Vw(t)||p < Ce{t~^ + 1) 
and the convergence statement for Vw^ follows. 



(61) 
□ 



Corollary 5.8 Consider fixed real numbers p E (|, 3) and a E (0, ^). Under the assump- 
tions of Theorem \5.3\ there exists a constant C depending onp,T, \\wo\\p, |||g|||o,p,T; llw^olli + 
||g||i.r and a, such that for all n E N, 



sup 7(t)-E(|K""(/i"'""-^)(t,x) -u(t,x)|) 
te[o,T] 



where ^(t) = {t'^p +t^ p^). 
Proof. For all {t,x) E [0, T] x M^, it holds that 
|K^"(/i'^'""'^)(t,x) -u(t,x)| 
< 



„ / (In n) 3 (In n) 3 1 

< C I V ^ + ^ r- + 



1 

n2' 



(Inn 



1 I 1 
9 / 



1 r— ^ 

K^" (/!"'--«) (t,x) - -^K,A^ - Xr-^"") A (x/j(<I>r'''))Mr,^5)l{.>a 

n 

- ^if,„(x - X,^'^"'^) A (Xi?,($r"'''))/^(T,X^)1{.>0 

(rr - y{t)) A x/?(</'(t))/i(^, x(O))P=-^((i0, dy, #) 
+ |K^"(w^")(t,x) - u(t,x)| 



(62) 



Ct 



with P^">^ as in Corollary 15. 5i By similar reasons as in (j56p . the expectation of the second 
term is now bounded by -^2M^^Rh. With the estimate in Theorem 15.21 we get that 

E |K,„ (/i"'^-«)(t, x) - u(t, :e)| <(L,„i? + M,J/i^C7(e„, ii, /i, T) 



+ 



2M,^ Rh 



n 



+ ||K^"(w^-)(i)-K(w)(t)||, 



Thus, from the estimates for and we deduce that for fixed p E (|, 3), 



£;|Ks„(/i"'^"'^)(t,j;) -u(t,x)| <C{1 + RhT){RhT) ^^"l'^'^ +CRh 



cln n) 3 



(cln n) : 



712 



/n 



+ ||w^"(t)-w(t)||H.i,P + ||K^"(w)(t) - K(w)(0||iy 
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where q = < We have used here again the Sobolev inclusions quoted in the proof 
of Theorem 13.71 and Lemma I3.3[ Now, by Lemmas 13.61 and 14.81 one has 

||VK^"(w)(t) - VK(w)(t)||,- <C\\ipe„ * w(t) - w{t)\\g < Ce„||Vw(t)||,- 



where we have also used part of Theorem 13.71 in the last inequality. On the other hand, 
||K^"(w)(t) - K(w)(t)||,- < Cy,„*w{t) - w(t)||p < Ce„||Vw(t)||p < Crhn, 

thanks to the estimate (|58p . From the previous estimates and Lemmas 14.91 and 15.71 we 
deduce that 

£1K..{,V«"''')«,:.0 -u(M-)l <C<5# + C<'°"'* 



'n 

__3_ 1 —l-Ll/i^i-i 

and the statement follows. [] 

6 Convergence rate under additional regularity assumptions 

Let us finally explain how the convergence rate can be slightly improved by assuming 
further regularity of the data wq and g. Since it is an adaptation of the developments in 
the previous sections, we only sketch the main arguments. 
First, it is possible to show that if the data wq and g are such that 

llttiollvy'n.p, sup ||g(t)||vF'"'P < oo, (63) 

te[Q,T] 

for some integer m > 1, then the mild solutions w'^,e > 0, given by Theorem 13.51 belong to 
the space Fm+i,p,r of functions v(f) such that 

m—l 



lll^'v|||o,p,T + |||^"^V||| 



i=l 

where -D* stands for the i-th order space derivative. To prove this, one easily first checks 
that wq belongs to that space, since the successive derivatives in the convolutions the 
heat kernel can be applied to the data wq and g. On the other hand, on can show by 
induction that the bilinear operators are continuous in F„_|_i^j,^t, and more generally, in 
the naturally generalized versions F^_|_i (^^.p) of the space Fi r^(T;p)- That is, the spaces of 
functions v such that 

m— 1 

ll|-C*v|||o,r,(T;p) + \lD"'^\ll,r,{T;p) 

1=1 

is finite. From this, one gets a local existence result in the space Fm+i,p,T, from which a 
regularity result can be obtained by arguments that can be adapted from those in the proof 
Theorem 3.2 in [10]. Moreover, one also checks that the norms |||w^|||m_|_i c^.p) are bounded 
uniformly in e > 0. 

Now, we impose additional conditions on the regularizing kernel (p, namely 
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i) /k3 ^{,x)dx = 1, 

ii) /jg3 \x\^+'^\ip{x)\dx < oo. 

iii) Jjg3 Xj^ . . . Xi^ip{x)dx = for all ii, . . . , v S {1, 2, 3} and r < m. 

Such function is called a cutoff function of order m + 1. Then, one has the following 
approximation result (see Lemma 4.4 in [23j): 



\iPs * w — w\\r < Ce 



m+1 ii^m+1 



WW 



for all w G W"^~^^'^ . Therefore, without any modification, for such function ip, the proofs 
of Lemmas 14.91 and 15.71 yield the same convergence results but at rate e™"*"^. 
By following exactly the same steps as in the previous section, we finally deduce 

Theorem 6.1 Assume the hypotheses of Theorems \5. 31 and moreover, that (j63p holds for 
some integer m > 1 and that ip is a cutoff of order m + 1. Then, we have for all i <n, 



E 



sup -Xi\ + - $j 



< C(r,u;o,g,a) 



1 1 

n2 (In nj 9 (Innj s 



and 



sup 7(t)£;(|K^"(/i"'^"'^)(t,x) -u(t,a;)|) 

te[0,T],x£M.3 



^ ( (Inn) 3 (Inn) 3 1 
< C I ^ , ' + ^ J. + 



n2 



m+1 

(lnn)~9~ 



where 7(t) was defined in Corollary \5. 81 where the constants now moreover depend on m. 

Acknowledgements: I would like to thank Mireille Bossy for suggesting to me the use of 
some techniques in [23j. 
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